Magnetic field influence on the proximity effect in semiconductor - superconductor 
hybrid structures and their thermal conductance 
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We show that a magnetic field can influnce the proximity effect in NS junctions via diamag- 
netic screening current flowing in the superconductor. Using ballistic quasi-one-dimensional (QID) 
electron channels as an example, we show that the supercurrent flow shifts the proximity-induced 
minigap in the excitation spectrum of a QID system from the Fermi level to higher quasiparticle 
energies. Thermal conductance of a QID channel (normalized by that of a normal QID ballistic 
system) is predicted to manifest such a spectral feature as a nonmonotonic behavior at temperatures 
corresponding to the energy of excitation into the gapful part of the spectrum. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 73.23.Ad. 
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The superconducting proximity effect is a mesoscopic 
scale phenomenon, which consists in the penetration and 
coherent propagation of Cooper pairs from a supercon- 
ductor (S) into a normal metal (N). The Cooper pair 
transfer into the normal metal can be equivalently de- 
scribed as an Andreev reflection process^ which con- 
sists of electron (with momentum p) conversion into a 
Fermi sea hole (with momentum — p) at the NS inter- 
face. The interference between an electron and the An- 
dreev reflected hole imposes a minigap onto the spec- 
trum of quasiparticle excitations near the Fermi level in 
the normal part of such a hybrid structure^, thus giving 
rise to pronounced features in its I(V) characteristics^!^ 
and thermoelectric properties^. Studies of the prox- 
imity effect have recently been made in various com- 
binations of materials, including junctions between su- 
perconductors and semiconductor structures^ support- 
ing a two-dimensional electron gas. In the case of elec- 
trons in a semiconductor structure weakly coupled to 
a superconductor, the minigap value discussed in the 
literature^sS, is much smaller than that of the 'mother' 
gap in the superconductor, both due to the mismatch 
vf vs between Fermi velocities in the two-dimensional 
gas [vp = {2Ep /my^'^] and the superconducting metal 
(vs), and also due to a possible Schottky barrier between 
them, with transparency 6 ^ g-2a/A (dependent on the 
length A of electron penetration into the barrier of the 
thickness a), Eg ^ ^^Ep < A. 

It has been noticed that the electron-hole interferences 
and the SN proximity effect survive at higher magnetic 
fields than the weak localisation - another quantum in- 
terference effecti^. This has been understood as a con- 
sequence of the fact that the interfering electron and 
Andreev-reflected hole retrace the same geometrical path 
in the normal metal, thus hardly encircling any magnetic 
fiux^. Therefore, another mechanism of magnetic field 
influence on the superconducting proximity needs to be 
taken into account, via a screening diamagnetic supercur- 
rent on the S-side of the hybrid structure. Since Andreev 
reflection takes place at the NS interface, where Cooper 
pairs flow, the incoming electron and the hole reflected 



by a moving condensate of Cooper pairs would be slightly 
shifted in momentum space; hence the ideal condition for 
them to retrace the same geometrical path is violated. 
As the orbital effect of the magnetic fleld on the normal 
metal or semiconductor side of the system is weak, the 
influence via diamagnetic screening may be the major 
factor of magnetic fleld influence on the superconducting 
proximity effect. 

Below, we analyze the influence of diamagnetic super- 
current in the system where the latter would be the only 
way a magnetic field might affect the proximity effect: a 
ballistic one-dimensional conductor connected in paral- 
lel to a superconducting bulk [Fig. la]. To be specific, 
we model such a conductor as a quasi-one-dimensional 
(QID) channel formed near the edge of a 2D electron 
gas in a heterostructure (x-y plane) with the side con- 
tact to a superconducting film, by depleting the 2D gas 
using a split top gate, and subjected to a weak magnetic 
field B = (0;0;i?). We show that the spectrum of low- 
energy quasiparticle excitiations in such a hybrid system 
has the minigap displaced with respect to the Fermi level 
to higher energies. 
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refiecting the fact that Cooper pairs in the channel are 
forced into the flow while tunneling from the bulk of the 
superconductor (where they are formed of two electrons 
with exactly opposite momenta) across the region of pen- 
etration of the magnetic fleld. [The Zeeman splitting ef- 
fect is also taken care of by the term aez (a is the spin 
projection) in Eq. Q]. As a result, each of the two elec- 
trons acquires the momentum shift 
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caused by the Lorentz force and equal to the differ- 
ence between the vector potential A = (0,A, 0) deep 
inside the superconductor, ^ = 0, and at its surface. 
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A = B5 tanh(L/2(5), where 5 and L stand for the London 
penetration depth and the superconductor fihn thickness, 
respectively. The spectrum described by Eq. iQJ can also 
be understood as one of the Bogolubov quasiparticlcs in 
the laboratory frame, where the equilibrium conditions 
are set by the heat reservoirs, for the condensate moving 
along the QID channel with the drift velocity Il/m (m is 
the effective electron in the semiconductor). According 
to Eq. 1^ the minigap is removed from the Fermi level 
when the field reaches the value 
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where Ed and ^ are the first critical field and the coher- 
ence length in the superconductor. 

The removal of a minigap from the Fermi level caused 
by a magnetic field would manifest itself in the trans- 
port properties of a hybrid system, such as the electron- 
mediated heat transfer. The ballistic quasiparticle spec- 
trum in Eq. gives rise to the thermal conductance 
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where kn(T) — irkgT/^h is the conductance of a nor- 
mal quantum ballistic wirei^i. At a zero magnetic field, 
the temperature dependence of k is activational, k(T < 
Eg/ks) oc e~^9/'"'^^, whereas at high fields, when there 
is no gap at the Fermi energy, k{T,B) — kn{T). The 
crossover from low to high fields takes place at B* 
[Eq. (j3Jl] and reflects the presence of a minigap Eg in the 
quasiparticle spectrum at finite excitation energies. This 
results in a nonmonotonous temperature and magnetic 
field dependence of the ratio k{T,B)/kn{T). 

The analysis of the quasiparticle spectrum formed due 
to multiple Andreev reflections in this paper is based on 
the standard weak-coupling approach to the proximity 
effect description in superconductor junctions with nor- 
mal metals and electron layers in semiconductors^. To 
be specific, we describe the QID confinement (provided 
by a gate) by the 2D electron wave function ip{x) lo- 
calized in the x direction, whose magnitude at the in- 
terface can be estimated from the boundary condition 
(p{0) — Xdx(p{0), with A standing for the electron pen- 
etration length into the barrier. The Fermi momentum 
of the QID system pp and 3D electron density on the 
semiconductor side are assumed to be much smaller than 
those in the superconductor, and we also take the tun- 
neling coefficient 9 ^ cxp(— 2a/A) as a small parameter. 
These assumptions enable us to neglect the influence of 
the normal system on the superconductor and to inves- 
tigate the proximity effect in the QID system without 
feedback. 

In the presence of a magnetic field B — (0; 0; B) it is 
convenient to choose the vector potential to be parallel to 
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FIG. 1: (a) Schematic view of a superconductor/QlD system 
junction, (b) Vector potential profile. 



the interface, A(x) = [0, A(x),0] in order to deal with a 
real order parameter in the superconductor. The vector 
potential A{x) acting on the normal electrons must be 
found self-consistently, taking into account the screening 
of the external magnetic field i? by a diamagnetic super- 
current^'^^. Inside the superconductor A(x) can be found 
from the London equation with the boundary conditions 
dxA{—a) — B and dxA[—a — L) = B as follows: 



A{x) = B6 



smh[{x + a + L/2)/S] 
cosh(i/2^) ' 



It is antisymmetric with respect to the middle of the su- 
perconductor: A{x = —a — L/2) — [Fig. 1(b)]. Since 
A(x) must be continuous at the surface of the supercon- 
ductor X — —a, in the semiconductor x > —a it varies 
as A{x) = B{x + a)+ BSta.nh{L/26). The width of the 
electronic wave function in the QID channel, Sx ~ kp^ 
and the barrier thickness a are both much less than 
L or 5, therefore, the vector potential acting on the 
QID electrons is virtually a constant: A{x) ~ A{—a) = 
BSta.nh{L/2S), which will be used below to determine 
the quasiparticle spectrum in the channel. 

We describe superconducting correlations in the QID 
channel using a pair of coupled equations for tppit) — 



annihilation operators: 
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where i? = 
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terizes the tunnehng coupHng to the superconductor 
and the electron momentum shift in the magnetic field 
n is related to the vector potential by Eq. ((SJ- In 
Eq. lO, G{t,t') = G{x = -a,x' = -a,t - t') and 
F(t,t') = F{x = —a,x' — —a,t — t') are the normal 
and anomalous Green functions of the superconductor 
at its boundary; a2 and as are Pauli matrices (tr* is 
transposed to a) . Since the size of the Fermi sea in the 
semiconductor wire is much smaller than in the super- 
conductor, one can ignore the dependence of G and F 
on the momentum parallel to the interface: only elec- 
trons in the superconductor moving nearly perpendic- 
ularly to the interface can tunnel into the QID wire. 
Since we are interested in the low-temperature regime 
fcsT ^ Eg ^ A, we will neglect the terms contain- 
ing the normal Green function G in Eqs. ((SJ). For the 
chosen gauge, the anomalous Green function of the su- 
perconductor, F in Eqs. Q has no phase factors, de- 
spite the presence of a magnetic field. For a weak field 
B <C Bel, its time Fourier transform can be estimated 
as F(e) « L-^Ep^ A/(A2 - + rf^J, with rj,,^ being 
the normal electron dispersion near the Fermi level in 
the superconductor. The integration over the perpendic- 
ular momentum gives F{e) « A/hvs{A'^ — e^Y^^ , thus 
giving us the minigap Eg — ■dF{e — 0) mentioned in the 
introduction and obtained in earlier publications^. 

The solution of Eqs. (O for e <C A is given by the 
Bogolubov transformation of the form 
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where bap and b[_ 



are Bogolubov's quasiparticle op- 
erators, and the excitation spectrum e„p is given by 
Eq. |lj (see Fig. 2). The Zeeman term in Eq. (P) 
turns out to be much smaller than the orbital one: 
Ez /vpli ~ g/kpT[\\iv{5,L) ^ 1 - unless the electron 
g-factor is anomalously large. 
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FIG. 2: Schematic view of the quasiparticle spectrum de- 
scribed by Eq. Q . 



Due to the motion of the QID condensate the exci- 
tation energy curve is tilted by energy Wi^IIsignp. The 
field B* [Eq. ©J at which the minigap is removed from 
the Fermi level is determined by the condition that 
VF^ = Eg. Note that at higher fields B* < B ^ B^, the 
quasiparticle spectrum remains gapful, with the center of 



the gap moved to energies 



Eg. 



Now we turn to the calculation of the thermal con- 
ductance k{T, B) of a long QID channel whose ends are 
kept at temperatures T and T + AT (AT < T). Since no 
heat can get into the strongly gaped superconductor, the 
middle of the wire represents a bottle-neck for the heat 
transport, so that we can analyze k{T, B) in the infinite 
wire geometry. The expression for the energy current op- 
erator jeijjt) in a wire can be found from the continuity 
equation dyjeijjt) — —dtPeiyt), where the density of en- 
ergy Peiyt) corresponding to the equations of motion Q 
is 
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where ipa{yt) = Ly Y.p'^a{t) eyip{ipy/h) with Ly be- 
ing the length of the QID system, andp = —iTidy. Using 
the Bogolubov transformation @ for V'q(^), for the den- 
sity of energy one finds 
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In order to satisfy the continuity equation with Pe{yt) 
given by Eq. © the energy current je (yt) must have the 
foUowing form: 
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In Eq. Q we have already omitted the terms containing 

bapb'La-p' and b^a-pbap' which vanish after the averag- 
ing. The averaged value of the energy current can be 
written as the sum of two contributions: 



-h j dp e+pdpe+p 
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The first of them can be attributed to the supercurrent 
flow and cannot transfer heat, whereas jq represents the 
heat current: 



= h- 



dp etpC'p'^tpnietp)- 
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The latter is determined by the energy distributions, 
n(e+p) and the group velocity, 9pe+p of quasiparticles. 
We express the energy currents (|1(JII and in terms 
of the "+" -branch of the spectrum ^ using the rela- 
tionship = — e^Q-p and the symmetry of the limits 
in the sum. The distribution functions of rightmovers 
{dpe^pQ > 0) and leftmovers {dpe'^p < 0) are assumed 
to be different and set by reservoirs, as n(e+p,T + AT) 
and n{e'^p, T), respectively. Using this, we determine the 
thermal conductance k{T, B) given by Eq. Q as the 
proportionality coefficient between the heat current and 
the temperature drop, jq = k{T,B)AT. 

Figure 3(a) shows the thermal conductance Q normal- 
ized by that of a normal wire as a function of ksT/Eg 
for different values of the magnetic field. Plot A is re- 
lated to i? = and shows how the conductance exponen- 
tially decreases at temperatures smaller than the mini- 
gap Eg. Curves B and C show what happens when the 




FIG. 3: (a) Temperature dependence of the thermal conduc- 
tance, K normalized by that of a normal wire, km for different 
values of magnetic field: (A) B/B* = 0.01, (B) B/B* = 0.95, 
(C) B/B* = 1.05, and (D) B/B* = 2. (b) Magnetic field 
dependence for different temperatures; (A) ksT/Eg — 0.1, 
(B) ksT/Eg = 0.3, and (C) ksT/Eg = 1. 



field crosses the value of B*, at which the edge of the 
minigap is about to reach the Fermi level. For B < B* 
(curve B), k(T)/ kn(T) is exponentially small only if 
ksT < Eg ~ i>i?n ^ Eg. When the temperature is in 
the interval Eg — vp^ < < Eg + wfII, quasipar- 

ticles with negative momenta p ~ —pp transfer heat, 
whereas the states with positive p are still unpopulated. 
This interval corresponds to the plato in curve B where 
the conductance k(T) is half of that in the normal state. 
At higher temperatures, ksT > Eg + vpll the asymme- 
try of the excitation spectrum no longer matters, and 
k(T) « kn{T). 



When the field exceeds B* (curves C and D), the de- 
pendence k{T)/ K]\r(T) becomes nonmonotonic. As in a 
normal wire, at low temperatures ksT <^ vp^ — Eg there 
are two left-moving and two right-moving modes capable 
of tranferring heat, which gives k{T) — K]\[(T). At inter- 
mediate temperatures vpH ~ Eg <^ ksT <^ Eg + wj^II, 
only the states with negative momenta contribute to 
the thermal conductance: k(T) = km{T)/2. At higher 
temperatures the conductance recovers a normal metal- 
lic behaviour. Finally, when B 3> B* the minimum in 
k{T) / KjsiiT) is less pronounced and the heat conductance 
behaviour becomes indistinguishable from that of a nor- 
mal wire. The magnetic field dependence of k/k^v is 
given in Fig. 3(b). 
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